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A new theoretical model for self dynamic response is developed using Vibration- Transit (V-T) 
theory, and is applied to liquid sodium at all wavevectors q from the hydrodynamic regime to the 
free particle limit. In this theory the zeroth-order Hamiltonian describes the vibrational motion 
in a single random valley harmonically extended to infinity. This Hamiltonian is tractable, is 
. . . evaluated a priori for monatomic liquids, and the same Hamiltonian (the same set of eigenvalues 

00 ' and eigenvectors) is used for equilibrium and nonequlibrium theory. Here, for the self intermediate 

' scattering function F''{q,t), we find the vibrational contribution is in near perfect agreement with 

, molecular dynamics (MD) through short and intermediate times, at all q. This is direct confirmation 

(N . that normal mode vibrational correlations are present in the motion of the liquid state. The primary 

transit efi'ect is diffusive motion of the vibrational equilibrium positions, as the liquid transits rapidly 
r-^ , among random valleys. This motion is modeled as a standard random walk, and the resulting 

theoretical F''{q,t) is in excellent agreement with MD results at all q and t. In the limit q — » oo, 
the theory automatically exhibits the correct approach to the free-particle limit. Also in the limit 
] 9^0, the hydrodynamic limit emerges as well. In contrast to the benchmark theories of generalized 

. hydrodynamics and mode coupling, the present theory is near a priori, while achieving modestly 

better accuracy. Therefore, in our view, it constitutes an improvement over the traditional theories. 

^ I PACS numbers: 05.20.Jj, 63.50.-fx, 61.20.Lc, 61.12.Bt 
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I. INTRODUCTION 



Important advances have been made in the theory of equilibrium thermodynamic properties of Hquids P, 0| • These 
' advances are characterized by the ability to calculate a priori the measured properties of real liquids, to an accuracy 
5 approaching the experimental accuracy itself [1, 0, IE] ■ The a priori nature of the theory is based on the key physical 
property of condensed matter, that the potential which governs the nuclear motion is given by electronic structure 
theory, in the form of the electronic groundstate energy as a function of the nuclear positions. In applying these broad 
theoretical foundations, two developments have been extremely helpful: (a) the development of pseudopotential 
' perturbation theory which gives effective internuclear potentials for nearly-free-electron metals in crystal and liquid 
phases @, and (b) the development of molecular dynamics (MD) computations to the point of providing highh^ 
[ accurate results, indeed capable of substituting for experimental data, when a good internuclear potential is used [Tj]. 
Based on these principles, a priori calculations have been made of binding energies of the elements [1] , thermodynamic 
properties of crystals [i| and liquids [T^ , and liquid static structure factors fTl|.ll2j| . What is important for the present 
work, MD has proven its reliability for no neq uilibrium properties, e.g. S{q,uj) for crystals and liquids [7., ,12.. 
^) • liquid resistivity [l^, and shear viscosity |l6l. [T7|. 
OO ', V-T theory was introduced when it was found that a theory based on two components of the atomic motion can 
' account for the equilibrium thermodynamic data of elemental liquids [T8| : (a) normal- mode vibrations in one (any) 
random valley, providing > 90% of the thermal energy and entropy, and (b) transits among a very large number of 
macroscopically equivalent random valleys, providing the remaining 10%. This theory is useful because the dominant 
vibrational motion is tractable, and its Hamiltonian can be evaluated for real liquids from electronic structure theory 
[T^ . To test this theory beyond equilibrium properties, it was applied to dynamic response in liquid sodium, where 
the vibrational contribution alone was found to give an excellent account of experimental results for the Brillouin peak 
dispersion curve [l^. A small correction for transits then produced agreement with MD data for the entire S{q,Lj) 



graphs 2]J. At this point we could see the possibility of a liquid theory for both equilibrium and nonequilibrium 
properties based on the same dominant (vibrational) component of the motion. Such a theory could prove useful 
because traditional nonequilibrium theories primarily describe the decay of fluctuations, through processes encoded 
in e.g. friction coefficients and memory functions, and these are concepts not present in equilibrium theories. Our 
point is not that the traditional description is wrong, but that a good part of it is already contained in the same 
vibrational motion that underlies the equilibrium theory. This view motivates the present application of V-T theory to 
self dynamic response. This application will provide a serious test of our theory, since self dynamic response has been 
thoroughly analyzed by the traditional theories of generalized hydrodynamics [l^l and mode coupling [2^. These 
analyses are discussed in liquid theory monographs, and have become a benchmark in dynamic response theories 

nil. 

The previous application of V-T theory to dynamic response was evaluated in the one- mode scattering approxima- 
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tion [2l|. Here, to work at larger q, it is necessary to use the full vibrational theory, including normal- mode scattering 
in all orders, i.e. keeping the displacement autocorrelation functions in the exponent. This formulation has not been 
studied previously, and what it reveals is remarkable to say the least. The vibrational contribution alone, in a priori 
numerical evaluation at all q, is in near perfect agreement with MD calculations through short and intermediate 
times, and the vibrational contribution alone also converges to the correct theoretical free-particle behavior at large 
q. Then, accounting for transit motion in leading approximation produces a theory in excellent agreement with MD 
calculations at all q and t. 

In Sec. IIA, the vibrational contribution to the self intermediate scattering function is derived, and its short-time 
behavior is examined, as well as its convergence to the free-particle limit. In Sec. IIB, the complete time dependence 
of the vibrational contribution is analyzed, and in Sec. IIC this contribution is compared with MD data. The transit 
induced correction to the vibrational contribution is derived and modeled in Sec. IID. The hydrodynamic limit is 
derived in Sec. HE. The complete theoretical results, vibrational plus transit, are compared with MD in Sec. IIIA and 
the two-step process by which theory approaches the free-particle limit is analyzed in Sec. IIIB. The salient theoretical 
features are summarized and discussed in Sees. IVA and IVB, and V-T theory is compared and contrasted with the 
classic benchmark theories in Sec. IVC. A brief sketch of the operational procedure of V-T theory may be found in 
the Appendix. 



II. THEORY 



A. The Vibrational Contribution 



We consider a system of N atoms in a cubic box at the density of the liquid, with periodic boundary conditions 
applied to the atomic motion. The atoms are labeled K = 1, . . . , TV and their positions are rj{{t) as functions of time 
t. The self component of the intermediate scattering function is [24, 25] 



K 



where the brackets indicate a motional average in an equilibrium state. The vibrational contribution expresses the 
motion in a single harmonic random valley extended to infinity [18] . In this motion each atom moves with displacement 
uxit) away from the fixed equilibrium position Rk, so that 

rK{t)^RK + UK{t). (2) 

Then the vibrational contribution to Eq. ([T]) becomes 



vib VI ' ^ / -p^ 



^ I vib 



where {■ ■ ■)yif, indicates an average over the vibrational motion in one (any) random valley. This is simplified by 
Bloch's theorem to 

F:^,{q,t) ^ l^e-2^-('^) e^-'-^-W 'J--(o)).... (4) 
The displacement autocorrelation functions are expressed in terms of the normal vibrational modes A (26j 

{q ■ UK{t) q ■ Mx(0)),„,b = 17 E'^^ ' "^^^^^ '^f^^!^^ (5) 

where T is the temperature, M is the atomic mass, wk\ is the Cartesian vector of the K component of eigenvector 
A, and uo\ is the corresponding frequency. The Debye- Waller factors are defined by 

WK{q) = \{{q-UK{Q)f)^^^, (6) 

and are given by Eq. ^ evaluated at t = 0. In addition to the vibrational averages, the right sides of Eqs. ([3]) and (21) 
are to be averaged over the allowed q- vectors at each g-magnitude, making Fmb{<l, t) a function only of q, as indicated. 
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Important properties of F^n,iq,t) are determined by its short-time expansion. To find this behavior we write, from 
Eqs. 0-®, 



kT , 1 

{q ■ UK{t) q ■ Mk(0))„, - 2WK{q) = — }2{q ■ w k x)^ — {cos u; xt - 1). (7) 

ivi — iU \ 

The expansion for ujxt << 1 is 

K 

The value at i = is 

F.Uq,0) = l, (9) 
which is the exact theoretical result, as can be seen from Eq. ([T]). The coefficient of t'^ is 

a{q) = kTq^/2M, (10) 
obtained with the aid of the eigenvector completeness relation [l^ 

^WKi,XWLj,X = SKL^ij, (11) 
A 

where i,j are the Cartesian directions. In studying the time dependence, it is advantageous to keep it in the exponent, 
as in Eq. ([8]), rather than to expand the exponential. The leading factor in Eq. ([8]) is the free-particle result, defined 
by 

F;_(q,^)=e-'^(«)*^ (12) 

Finally, the coefficients of in Eq. ([5]) are 

bK{q)^'^Y.(^-^Kxru;l (13) 

where the unit vector q has been introduced so as to factor out q^. Without approximation, the i^-dependence of the 
bxiq) cannot be ignored. 

In dynamic response, free-particle motion becomes dominant at short times and distances. Free-particle behavior 
is contained in the self dynamic response [U, [2^, and should emerge in the large-g limit of Eq. ([4]). Since the right 
side of Eq. ([7]) is negative and proportional to , as q increases -FJj^ (g, t) will drop off more rapidly with increasing 
t. Therefore at high q only the lowest-order term in Eq. ^ will be relevant, so that 

lim F^^,{q,t) = F^„^^iq,t). (14) 

q — too ^ 

The process of the approach to this limit will be examined below. 

B. Free, Intermediate, and Convergence Periods 

We calculated F^^^{q,t) for the 17 q values listed in Table I. Representative graphs are shown in Fig. 1. Each 
curve has the same characteristic shape; it starts at 1, decreases on a uniform (g-indcpendcnt) timescale, then levels 
off and converges to its long-time limit. The displacement autocorrelation functions, Eq. ([5]), decay to zero as time 
increases. This decay is called the "natural" decorrelation, and is responsible for the entire time dependence of 
Fyibiq^t)- Analysis reveals three distinct periods in the curves. 

In the initial period, the atomic motion is free. During this period, F^^f^{q,t) is dominated by the leading factor in 

Eq. ([HI, e~°^'^*^. To estimate the duration of this period, let us average the right side of Eq. over the atoms, and 
over the star of q; then with the eigenvector orthonormality relation 

^WKi,XWKi,X' ^ Sxx' (15) 
Ki 
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1 2 3 

time (ps) 

FIG. 1: Fyitiq, i) for a wide range of q, showing the initial decrease on a g-independent timescale, and universal behavior where 
the function converges to its t = oo limit. 



we find 



Accordingly 



1 

1 3N' 



(16) 



(17) 



where 



kTq^ 
24M 



(18) 



and where (i^x) = {•^N) ^J2x^x- ^^^^ motion period will last until the term in Eq. P7|) begins to be felt; 
let us therefore choose Tf, the duration of this period, as the time when b{q)t'^ = 0.1a(g)t^, giving r/ = 1/1.2/ {mj). 
Because of the decoupling approximation used in deriving Eq. (|17p . Tf is independent of q. {i^x) is related to the 
characteristic temperature 6*2 by | (^{hwx)'^^ = (^6*2)^ 1^9]. For our liquid Na system 62 = 154.1 K [23|, and we find 
Tf = 0.070 ps. Our calculations show that F^ii,{q,t) begins to depart from e""'-*^*^ at a time near Tf, specifically at 
around 0.08 ps at g = 0.30 a^^, and decreasing to around 0.05 ps at g = 3.50 Gq^. 
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TABLE I: Infinite time limit of F^ifj{q,t) and decay factors 7(g) of the transit induced decorrelation, according to Eq. (1281) . for 
different values of the wave vector q. 









.29711 


0.94185 


0.1733 


.70726 


0.71699 


0.9222 


.91575 


0.57671 


1.4676 


1.0148 


0.51138 


1.7500 


1.0917 


0.46182 


1.9753 


1.1050 


0.45332 


2.0146 


1.1443 


0.42926 


2.1312 


1.2547 


0.36482 


2.4595 


1.5052 


0.24047 


3.1805 


1.7577 


0.14876 


3.8169 


2.0041 


0.08776 


4.2975 


2.2529 


0.04903 


4.6103 


2.5064 


0.02573 


4.7423 


2.8498 


0.00985 


4.6527 


3.2000 


0.00339 


4.3396 


3.5008 


0.00126 


3.9994 


6.0013 


0.00000 


4.2271 



The intermediate period in Fig. 1 follows the free period. Here the strong decrease of F^^j^{q,t) continues, but the 
function is not approximated by the e""'-*^* factor in Eq. ([8]). In terms of the power series expansion of the right 
side of Eq. ([7]), increasingly higher orders contribute while retains a smooth t-dependence. This property is 

q- independent. 

At some point, the set of cos Loxt dephase and begin to cancel, starting from the highest frequency and continuing 
to the lowest, which is the last to dephase. Since the highest frequency in our system is ujmax — 25.5 ps~^ [2l|, the 
dephasing will begin around 27: /uJmax = 0.25 ps. From Fig. 1, this is close to the uniform (q-independent) time when 
the initial decrease ends and begins to converge to its long-time limit. Hence we associate the convergence 

period with the normal-mode dephasing process. 

The t — > 00 limit is obtained by setting to zero the displacement autocorrelation functions in Eq. ^ : 

^J..('Z,-) = ^E^"'^"^'^^- (19) 

K 

In approaching this limit, when the displacement autocorrelation functions are sufficiently small, a first-order expansion 
of the time dependent part is useful. With Eq. ([5]) this expansion reads 



F^^,{q,t) = F^^,{q,^) 

M ^-^ Uu. 

K \ 



From Eq. (I19|) . it is the Debye- Waller factors which set the level to which F^^^^{q,t) decreases in Fig. 1, and after that 
decrease, Eq. (PH)) applies. Values of F^-f^{q, 00) are hsted in Table I. 

Fig. 1 shows a remarkable similarity of the curves in the convergence period. The main feature is a broad minimum 
around 1 ps, where F^-f^{q, t) lies below F^.^{q, 00), and a final increase of F^-j^^q, t) to arrive at F^^^{q, 00) by around 
2 ps. These timings are accurately independent of q. Furthermore, superimposed on this broad shape is a set of 
small oscillatory features whose timings are also accurately independent of q. This property can be understood with 
the aid of a small approximation in Eq. (|20[) . In a single random valley, the atomic sites are all inequivalent, for the 
same reason that different sites in a crystal unit cell are inequivalent. If we neglect the coupling of this inequivalence 
between the eigenvectors and Debye- Waller factors, we can average the (q • wkxY as in Eq. (fT6|) and transform 
Eq. (HOD to 



kTq^ 1 cosLU\t 



(21) 



The function in brackets now exhibits uncoupled dependence on the variables kTq^/M and t, where the t-dependence 
is parameterized by the set {uj\} of normal mode frequencies. Numerical tests verify that Eq. (|2ip is rather accurate 
in the convergence period, and this explains the uniform ((/-independent) time dependence of the curves in Fig. 1. 







0.1 0.2 
time (ps) 



0.3 



0.4 



FIG. 2: For F''{q,t) the vibrational contribution (lines) and the MD data (dots) agree extremely well to around 0.2 ps (fine 
vertical line). 

The properties of F^^f^{q, t) were also examined for several different random valleys in our system, and only insignifi- 
cant differences were found. Hence the uniformity of random valleys previously found for equilibrium thermodynamic 
functions [23 is extended to the time correlation function studied here. 



C. First Comparison with MD 

Fig. 2 shows superimposed graphs of our MD results, Fljjj{q,t), and the vibrational contribution, F^^i^{q,t). At 
each q, the MD and vibrational graphs are in near perfect agreement up to an "initial departure" at around 0.2 ps. 
Beyond the initial departure, F^ii,{q,t) turns away from F^j]j{q^t) and proceeds to converge to its long time limit. 
Fig. 2 shows a remarkable constancy of the initial departure for all q, and even a remarkable constancy in the shapes 
of the curves in the vicinity of the departure. 

Let us consider what Fig. 2 means for our analysis. F^^f^{q, t) is based on pure vibrational motion in a single random 
valley. On the other hand, Ff^j]j{q,t) is based on the real liquid motion, which has both vibrational and transit 
contributions. From Fig. 2, the vibrational motion completely dominates f^£,(g,i) up to the initial departure. After 
that, the difference between F^^^(g,t) and F^-f^{q,t) is due to transits, which are present in the MD data but not 
in the vibrational calculation. However, even though transits are going on continuosly in the MD system, at a very 
high rate because the temperature is around T„n the effect of transits does not appear immediately in F^jj[q,t)^ 
but only after the motion departs from the pure vibrational motion. This observation will be used to calibrate our 
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decorrelatioii model in the next Section. 



D. Transit-Induced Decorrelation 



Our goal here is to model the effects of transits in a way that is consistent with the a priori determined vibrational 
motion in a single random valley. We start by returning to the definition, Eq. ((T|), which becomes 

F^^(g,i) = l//^e-'^-(-W--(°))\ \ . (22) 

\ \ I vib' trans 

This form expresses the fundamental insight of V-T theory that the motion of the atoms consists of vibrations that 
are periodically modified by transits, allowing us to consider the effects of the two kinds of motion as two consecutive 
averages. Expanding each tk and evaluating the vibrational average as in Eqs. ©-([l]), we find 

Fyj.{q,t) = i- /e"''J-(^^^(*)"-^^('^)) g-2WK(q)g(q-Wif(t) q-UK-(0)>„,A _ ^23) 

^ — ^ \ / trans 

K 

Notice that the equilibrium positions still carry time dependence because we have not yet evaluated the transit average. 
From this point of view we identify two ways in which transits will modify F''{q, t). The first arises from transit-induced 
changes in the atomic equilibrium positions Ridt), while the second way is through their effect on the displacement 
autocorrelation functions {q ■ uxit) Q ■ 'f'KiO))^^,- Given the fact that transits change the equilibrium positions on 
very short time scales compared to the vibrational motion, we expect the first effect to be largely decoupled from the 
second, so we make the approximation of full decoupling and find 

^) ^ i_ /g-*q-(RK(t)-H/c(0))\ /^-'2WK{q) {q-UKit)q-UK{Q)),,A (24) 

trans \ / trans 

Now we can consider each transit average separately. 

Let us abbreviate the first transit average as Ax{t). In V-T theory the motion of RK{t) is entirely responsible for 
self diffusion, hence it is appropriately modeled as a random walk. The single atom transit rate v is the number of 
transits per unit time in which a given atom is involved. Consider an increment 5t sufficiently small that an atom is 
very unlikely to be involved in more than one transit. Then in 5t, Axit) changes by 



MK(t) = ( e-*''-«^^*+''*^ -e-^'J-"^^'^ e*«-^^(°)\ . (25) 

\ L J / trans 

In 5t, each atom transits once with probability v5t, or else does not transit. If atom K does transit, Rxit + St) = 
Rnit) + 5Rk- If atom K does not transit, Rxit + St) = Rxit). Eq. (|25p becomes 

SAK{t) = /[e-"J-*-f^^- - 1] e-'qiRKit)-RKm\ ^st. (26) 

\ / trans 

We assume |(5i?i<-| ~ SR, the same for all transits, while the direction of SRk is uniformly distributed and uncorrelated 
with the other factors inside the sum. Then — 1] can be separately averaged over angles and Eq. (|26|) can 

be written 



SAK{t) 



St 

where 

sin qSR 



j{q)/e-'^i-iR'<W-R^m\ (27) 

\ / trans 



(28) 



qSR 

Since the transit average on the right hand side of Eq. ([?7|) is AK{t) and Ak{0) = 1, the equation integrates to 

AK{t) = e-T^-?)*. (29) 
Finally, from this result and Eq. (g]) for F^^f^{q,t), Eq. (pi)) becomes 

F^^(g,i)=e-^(«)*(FJ,,(g,i)),_,. (30) 
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time (ps) 
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FIG. 3; Comparison of FYj,[q,t) (lines) with Fl4]j{q,t) (circles) for various q. The crossover bump near Tc is enlarged in the 
graph for q — 1.51 a^^ . 



So the random walk for Rxit) in Eq. yields the last relation, with the displacement autocorrelation functions 
in F^^fj{q,t) still having arbitrary time dependence. 

Now we must account for the fact that the transit-induced decorrelation does not begin to operate until some 
time has passed (Fig. 2). We do this with the simplest possible model, replacing the decay factor e"^'*-** by the 
decorrelation function D{q,t), defined by 
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0.02 



0.00 

5^ 



-0.02 



1 2 3 4 5 6 7 
time (ps) 

FIG. 4: Deviations of Fvj'{q,t) from Flj]j{q,t), Eq. (|35|l . The dotted line is q — 0.30 ^, the four dashed hnes are (from the 
lowest) q = 0.71, 0.92, 1.01, andl.25 a^^ , and the solid lines show negligible deviation for q > 1.51 a^^ . 

where is the crossover time, to be calibrated. Then from Eq. (|30p . om' transit decorrelation model is now 

F{.^{q,t) = {F:^,{q.t)),^^^^D{q,t). (32) 

In the random walk model, the self diffusion coefficient is Z? = ^^{SR)^ . A measurement of 5R was obtained from 
MD simulations at 30 K [2^, and the result can be used here because 5R has only weak T dependence. On the other 
hand, v \s & strong function of T, so we use the MD evaluation Dmd = 5.61(10~^cm^/s) at the temperature of the 
present study, 395 K [2§|. For comparison, the experimental value for hquid Na at 395 K is 5.08(10~^cm^/s) [soj . 
The results for 6R and ly are 

SR = 1.75 flo, 
ly = 3.9 ps"^ 

The value of v is not far from our previous estimate of 2.5 ps~^ [3l|. Data for 7(g) are listed in Table I. 

The second effect of transits is felt on the displacement correlation functions. As noted in Sec. IIB, these functions 
already decay to zero by the natural decorrelation among normal modes. As a model, we might suppose that transits 
cause sudden vibrational phase shifts, and therefore enhance the decay of the displacement autocorrelation functions. 
This enhanced decay will cause F^^i^{q,t) to decay more quickly to F^jj((7, 00). Referring to the curves in Fig. 1, the 
main effect will be to slightly raise the broad minimum around 1 ps, and bring the curve to F^j^[q, 00) noticeably before 
2 ps. This is an interesting physical effect, but is small; in fact it is extremely small compared to the decorrelation 
caused by the motion of the equilibrium positions, as modeled in Eq. (|3T|) . We therefore neglect the transit-induced 
decorrelation of the displacement autocorrelation functions in the present study, with the final result that 

F^T{q,t) ^ F^,,{q,t)D{q,t). (33) 

Now the only thing that remains to be determined is the parameter Tj,. The foremost property revealed in Fig. 2 is 
that Tc must be independent of q. For this reason it is conceivable that vibrational information alone can determine 
Tc- We have clues to this effect - for example, Tc is close to the start of dephasing at 0.25 ps (Sec. IIB) - but we have 
no certainty. We therefore calibrate Tc by comparison of the theoretical Fyrp^q^t), Eq. ([55]) . with the MD results in 
Fig. 2, at t around Tc- This gives Tc = 0.22 ps. The reason Tc is a little larger than the initial departure time of 0.20 ps 
in Fig. 2 is to allow the discontinuous D{q,t), Eq. ([3T|) . to better fit the actual smooth crossover behavior. When Tc 
is determined in the same way for different random valleys, the scatter in Tc is around ±0.01 ps. 

E. The Hydrodynamic Limit 

The hydrodynamic limit corresponds to g ^ and t larger than a characteristic time. As q 0, from Eqs. ([4|-([6|. 
Fvibi^^t) 1 + 0{q^) for < t < cx). Also as g ^ 0, 7(5) 0{q^), and from Eq. D{q,t) e-^(«)*[l + 0{q^)] 
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FIG. 5: Comparison of Svt{<1j'^) (lines) with Sliu{q,uj) (circles) for various q. 



for t > Tc. Then from Eq. dSS]), F(,j,{q,t) is e"''(«)*[l + 0(g^)]. We drop the term in 0{q'^) and take the g ^ limit 
of 7(g), Eq. (I28|), to find 



lim F^g.(g, t) = e"-""?'* for t > 



(34) 



where we used the random walk model expression for the self diffusion coefficient. The right side is the correct 
hydrodynamic limit of F''(q,t). 
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FIG. 6: Comparison of theory (solid line) and MD (circles) for the normalized halfwidth of S''{q,u!) for all q listed in Table I. 
The dashed line is the free-particle limit. 



III. COMPARISON OF THEORY WITH MD CALCULATIONS 



A. Self Dynamic Response 




FIG. 7: Comparison of theory (solid line) and MD (circles) for the normalized value of S"{q,0) for all q listed in Table I. The 
dashed line is the free-particle limit. 

Graphs comparing theory with MD arc shown in Fig. 3 for q = 0.30 — 2.51 a^^ ■ The time required to decay to 
zero decreases strongly as q increases, from around 28 ps at g = 0.30 a^^ to 0.6 ps at g = 2.51 Uq^. The decay rates 
7(g), Table I, show a corresponding strong increase as q increases. Since the decorrelation begins at Tc = 0.22 ps, 
nearly the entire curve at q = 0.30 Qq^ is the decay process. The decorrelation steadily becomes a smaller part of the 
curve as q increases, until it affects only the tail of the curve at g = 2.51 Og ^. These changes are accompanied by a 
change in shape of the graphs, from a near-exponential at q = 0.30 a^^, going over to a qualitatively Gaussian shape 
at q = 2.51 ag ^. All of this behavior has the important consequence that the decay rate ^{q) eventually becomes 
irrelevant as q increases. 

Generally speaking, the crossover bump in the theoretical curve around Tc cannot be seen in graphs as small as 
those in Fig. 3. The prominence of the bump depends on both the slope and magnitude of the curve at Tg. It is most 
prominent at g = 1.51 a^^, where it is enlarged in the inset. The deviations 



AF%q,t) = F^T{Q,t) - Fljij{q,t) 



(35) 



12 



for representative q are shown in Fig. 4. The largest deviations in our study are included in the figure. For most 
q and most t, the deviation is < 0.005 in magnitude, which we consider insignificant. Larger deviations are seen 
in two places, (a) The positive spike at t « 0.22 ps is due to the crossover model, Eq. ([3T|) . This error could be 
removed by smoothing, but since it is so small we are willing to forego the introduction of a smoothing model, (b) 
The only significant deviations arc those in the negative dip at t ~ 0.3 — 2.0 ps. These deviations reach magnitude 
~ 0.01 — 0.02 for the dotted curve {q = 0.30 a^^) and the four dashed curves {q = 0.71 — 1.14 Qq^). Notice the 
shape at q — 0.30 Oq ^ is the same as the others, but stretched out over a much longer time. This error arises from 
the broad shallow minimum in F^^f^{q,t) observed in Fig. 1. That minimum is not entirely in error, however, since a 
weakened copy of it is present in F^j-, (g, t) . From a detailed examination of the data, we conclude that the error can 
be eliminated by adding a theoretical model for the transit-induced decorrelation of the displacement autocorrelation 
functions. This decorrelation process was discussed at the end of Sec. IID. 

Graphs comparing theory with MD for S'^{q,Lo) are shown in Fig. 5. The only significant error corresponds to the 
negative dip in AF^{q,t) for a limited q range, in Fig. 4. This q range is near the location of the first peak in S{q), 
at qm = 1.05 Og ^ in our system. 

Further comparison of theory and MD is shown in Figs. 6 and 7. The quantity graphed in Fig. 6 is Auj{q)/ Dq'^ , 
where Auj{q) is the halfwidth of S''{q,uj) and Dq^ is the halfwidth in the diffusion limit 0, H^. Fig. 7 shows the 
quantity TrDq'^S''{q,0), where {nDq'^)^^ is S"*(q, 0) in the diffusion limit In these graphs we used the same value 
Dmd for Na at 395 K as in calibrating 7(5) following Eq. ([55)1 . These figures will allow a comparison of the present 
theory with the traditional generalized hydrodynamics and mode coupling theories in Sec. IV. Here we observe that 
Figs. 6 and 7 show a rather good agreement between theory and MD for the entire q range, including the approach 
to the free-particle limit. Further, most of the error revealed in Figs. 6 and 7 appears in the vicinity of q^, and is 
attributed to our neglect of the transit-induced decorrelation of the displacement autocorrelation functions in Eq. (4) . 



According to Eq. (fH)) . F^^^^{q,t) reaches the free-particle limit as q ^ 00. When the transit-induced motion of the 
equilibrium positions is accounted for, the V-T theory expression is Eq. (l33t . with D{q,t) given by Eq. (jST]) . The 
approach of Fyj,{q, t) to the free-particle limit is therefore to be understood by watching the function F^^i^{q, t)D{q, t) 
as q increases. The process contains two steps, both operating at all times, but in effect more or less sequential, as 
follows. 

(a) In the first step, because F^^^{q, t) decreases to zero in an ever shorter time as q increases, and because D{q, t) — 1 
for t < Tc, the decorrelation function effectively approaches 1 as q increases. This process is apparent in Fig. 2. Let 
us define qc by the condition that D{q,t) can be replaced by 1 for q > qc. Then 



B. Approach to the Free-Particle Limit 
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FIG. 8: F''{q,t) at q ~ qc. Solid line is V-T theory, circles are MD, and dashed line is the free-particle theory, Eq. (|12|l . 
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(b) In exact free-particle motion, the self intermediate scattering function is given by Eq. p2|) . But when step (a) 
is completed and Eq. ([M)) holds, and F^it,{q,t) agrees with Fljj-,{q,t), these functions have not yet reached the free- 
particle limit. This is also shown in Fig. 8, where the theory and MD curves differ noticeably from the free-particle 
curve. The situation exemplifies the point discussed in Sec. TIB, that in the intermediate time period, F^ii,iq,t) 
continues its steep decrease but is not well approximated by the e""^"^^* factor. So the final step in arriving at the 
free-particle limit is to increase q beyond qc, until F^^i^{q,t) concides with Fj^^^{q,t). 

Fig. 9 shows the Fourier transform of Fig. 8. Again the theory and MD are in agreement while both differ from the 
free-particle curve. Another view of the approach to the free-particle limit is seen in Fig. 6. For the halfwidth, the 
difference between theory and MD is ar ound 0.8% for q = 2.51 - 3.50 a^^ , and is 0.1% at g = 6.00 a^^ . This level 
of agreement is in strong contrast to the much larger difference of theory and MD from the free-particle halfwidth. 
The same property for S^{q, 0) is observed in Fig. 7. Clearly theory and MD are in substantial agreement long before 
they arrive at the free-particle limit. Then, the final approach to the free particle limit is determined entirely by the 
vibrational contribution, and is very slow, possibly algebraic, as seen in Figs. 6 and 7. 

IV. SUMMARY AND CONCLUSIONS 

A. Vibrational Contribution 

The main conclusions of this study are expressed through a discussion of the self intermediate scattering function. 
We begin by comparing properties of Fyu,{q,t) with exact theory and with MD results. 

(a) The short-time expansion is Eqs. (H]), pO)) . and p3|) . The value at i = is F^^^{q, 0) = 1, the exact result. The 
expansion gives the free-particle limit F^^^^{q,t) = e^"'^'''* as the leading factor in F^-f^{q,t), at all q. The result is 
not trivial. It appears because the eigenvector completeness relation Eq. ([TT]) decouples the free-particle motion in 
the exponent. This property ensures that V-T theory will automatically give the free-particle behavior as q — > oo. 

(b) At the end of the period of free-particle motion, at t = t/, F^^i^{q,t) starts to depart from e^^^'''* . F^^ij{q,t) 
continues its strong decrease throughout the intermediate period, until F^-f^{q,t) levels off and starts to converge to 
F^jj((7, oo). The power series in does not usefully represent F^^^{q,t) in the intermediate period. F^^i^{q,t) is in 
very good agreement with Fljjj{q,t) for all times up to near the end of the intermediate period, and all q, Fig. 2. 
This period of agreement between F^^f^{q,t) and F^jj-,{q^t) is three times longer than the duration Tf of free-particle 
motion, providing direct confirmation that normal mode vibrational correlations are present in the motion of the 
liquid state. 




FIG. 9: S''{q,uj) at g ~ Qc. Solid line is V-T theory, circles are MD, and dashed line is the free-particle theory. 
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B. Transit Contribution 



To complete the theory for the self intermediate scattering function, we must include the transit motion. What 
is needed is a physically motivated model which will not interfere with the above listed accurate properties of the 
vibrational contribution. In Eq. ^ for F^^i^{q,t), there are two places where transits will cause decorrelation of the 
purely vibrational atomic motion. 

(a) The first effect of transits is to make the equilibrium positions time-dependent in Eq. i.e. Rx{t). This 
motion is modeled as a random walk of transit jumps 5Rk, all of the same magnitude 5R, and uniformly distributed 
over angles for each atom. This transit motion then decouples from the vibrational motion, to give the complete 
theory in the form of Eqs. (|3ip and (|33p . The random walk model is calibrated a priori from MD data for the transit 
jump distance and the self diffusion coefficient. The crossover time Tc is the only parameter calibrated with the aid of 
MD data for self intermediate scattering. The transit decorrelation embodied in D{q, t) constitutes a major correction 
to F^-f^{q,t), and brings Fyrp{q,t) into excellent overall agreement with F^^^(g,t). 

(b) The displacement autocorrelation functions, written in Eq. ([5]), decay to zero as t — > cxd, by virtue of the 
natural decorrelation among normal mode vibrations. The second transit effect will be to enhance this decay. To 
make a proper model for this will require additional study, especially regarding the form of the decorrelation and the 
importance of the K dependence. In the meantime, the present study shows that this transit effect is small enough 
to neglect entirely and still have a highly accurate theory. 

For the self function, the transition to free-particle behavior involves two sequential steps. First, with increasing q, 
F^^^{q,t) goes to zero in a shorter and shorter time, until at qc F^^^{q,t) vanishes for t > Tc. Then for q > qc, D{q,t) 



in Eq. 



is effectively 1, so that Fyrp{q,t) = F^^f^{q,t), and Fyrp[q^t) agrees with Fljjj{q,t). As q increases from 



qc, theory and MD remain in agreement, and together they approach the free-particle limit, 
within the vibrational contribution. 

In contrast, the hydrodynamic limit depends entirely on the transit decorrelation. As g — > 
for t > Tc, and this produces the hydrodynamic limit in Eq. 



This final approach is 
0, F{.^iq,t) - 



0.15 




FIG. 10: Error in di.ui/Dq^ as function of q/qm- Line is present theory, {qm = 1.05 ay ^), open and filled circles are Ar from 
[2^ and [2^ ] respectively, and filled diamonds are Rb from ^3]. 



C. Comparison with Traditional Theories 

The self dynamic structure factor for liquid Ar near the triple point was studied by Levesque and Verlet [2^ , who 
compared several generalized hydrodynamics models with MD for the quantities graphed in Figs. 6 and 7. A similar 
study was made for liquids Ar and Rb by Wahnstrom and Sjogren who carried out their analysis in mode coupling 
theory. These studies are discussed by Hansen and McDonald [23| (p. 266) and by Balucani and Zoppi [25] (Sec. 5.3), 
and have become the benchmark theories for S'^{q,uj). Levesque and Verlet use the memory function formalism "to 
give a simple phenomenological fit for the computed self intermediate scattering function." The fitting process uses 
the q-dependent coefficients of t^,^^, and in F''(q,t), and is also adjusted to give the correct ideal gas (high-q) 
limit. Wahnstrom and Sjogren split the memory function into a rapid binary collision part and a slowly varying 
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FIG. 11: Error in nDq^ {q,0) as function of q/qm. Line is present theory, (qm ~ 1.05 a^^), open and filled circles are Ar 
from [52] and [2^ respectively, and filled diamonds are Rb from [2^ . 



collective recollision part. Their calibration also uses the coefficients of t'^,t'^, and in F^(g, t), and additional longer 
time information. The application of V-T theory constitutes an extreme contrast to these traditional theories. The 
collision concept does not appear in V-T theory. The vibrational contribution alone, a priori and for all q, accurately 
accounts for the short- and intermediate-time behavior of F''(q,t) and accurately accounts for the free-particle limit 
as well. Hence F^^^{q,t) provides much of the information used for calibration by the traditional theories. Then, only 
the transit contribution remains to be addressed. 

Since we are proposing a radically new theory of self dynamic response, we should compare numerical accuracy with 
the traditional theories. The standard for this comparison is the error of theory from MD for the quantities graphed 
in Figs. 6 and 7, the width and height of S^^q.uj) [2I [H, [H ill . These errors, in the form of (theory-MD)/MD, are 
compared for traditional and present theories in Figs. 10 and 11. The Levesque and Verlet data for Ar are from their 
tables [121, and the Wahnstrom and Sjogren data for Ar and Rb are read from their graphs [1^. The large error for 
Rb at large q is attributed to a difference of the self diffusion coefficient between theory and MD [23] . In the overall 
comparison, the error in the present theory is smaller than in the traditional theories. However, in our view, more 
important than the comparison of numerical accuracies is the near a priori character of the present theory. Beyond 
the standard diffusional model for the motion of the vibrational equilibrium positions, only the single scalar parameter 
Tc is needed to calibrate Fyrp{q,t) for all q and t. 



APPENDIX A: OPERATIONAL PROCEDURES OF V-T THEORY 



Since the operational procedures of V-T theory are somewhat novel, a brief summary will perhaps be useful. One 
begins with an interatomic potential for the system of interest. An MD system is constructed, and quenches are made 
from the liquid to minimum potential energy structures, where the equilibrium positions and the dynamical matrix 
are evaluated. Except for three translational eigenvalues, which are zero to numerical accuracy, all normal mode 
eigenvalues must be positive. In the first application to a given system, one needs to make some tests to identify the 
random structures (as differentiated from symmetric structures) , and to verify that the random structures do indeed 
dominate the potential energy surface [27l. [29l [33 | . This step is to check the "single random valley" approximation, 
whose verification is still in progress [27[. From this point one can proceed by working with a single random valley. 
The dynamical matrix is diagonalized to find the frequencies ux and eigenvectors wkx- The vibrational contribution 
to any thermodynamic function [isl . [l9| or to a time correlation function [26] can be expressed in terms of these 
quantities. This is illustrated by the equations of Sec. II. Numerical evaluation of sums /-^ proceeds by direct 
summation over the iN — 3 modes with nonzero ujx- At the present stage in the theoretical development, the transit 
contribution to a statistical mechanical average is accounted for by a macroscopic model, such as the one reported 
here in Sec. IID. 
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